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We show that at a second order phase transition, of (f) 4 like system, a necessary 
condition for stretched exponential decay of the time dependent structure factor is 
obeyed. Using the ideas presented in this proof a crude estimate of the decay of the 
structure factor is obtained and shown to yield stretched exponential decay under very 
reasonable conditions. 
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The decay of disturbances in some physical systems is known to be a stretched 
exponential in time. The canonical list of those systems includes glasses, [1-3] 
polymer solutions [4,5] and such a form of decay is also characteristic of dielectric 
and viscoelastic relaxation. [6,7] In recent work [8,9] we have shown that this form 
of decay is much wider spread than expected. We claim, in fact, that such a decay is 
present in many non-linear systems such as magnets at their critical point, systems of 
the KPZ type etc. These systems share the property that a disturbance of wave vector 
q , decays at equilibrium with a characteristic decay rate, 0) q , that is proportional to 

g M with ji > 1 . Since, the problems involved are very difficult, the derivation has to 
rely on approximations. Therefore, it is very important to have some exact results that 
may also suggest an intuitive understanding of the stretched exponential decay. We 
start by giving an exact proof that the condition for a slow form of decay is fulfilled at 
the critical point of a system of the (j) 4 type. We continue then, by employing the 
same basic ideas, to argue that indeed the decay is a stretched exponential in time. 
The systems we consider here are very different than those considered so far in the 
context of stretched exponential decay and we are interested in q dependent 

correlation functions, yet, the final step in the derivation is similar in its mathematical 
structure to other derivations in which stretched exponential decay has been obtained. 
[10-13] 

The system we have in mind is a classical (j) 4 system at its critical point or any 
generalization of it. We will also assume that the classical Hamiltonian, W , that 
gouverns the statics is even in the field . It is well known that Langevin dynamics 
of the system leads to a Schroedinger like evolution equation [14-17] 
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at 

where y/ is related to the probability distribution of obtaining a given field 
configuration {<p q } at time t, P{<t> q ,t} , by P{(/> q ,t} = exp[-W /2ty{(f) q ,t} and where 
the "Hamiltonian" 

II - -> XA + } (X =X„>0) (2) 

^ qX d^_ q 2d^_ q 4d<j> q d<j>J " " "* 

is Hermitian, non negative and its only eigenfunction with eigenvalue zero, is the 

y 

ground state y/ c = Pe q n °cexp[-W/2] . In the following we will chose X q to be 
proportional to q e with 6 that is not too negative. (The meaning of the last phrase 
will become clear later). Past studies of the dynamics of the (j) 4 system concentrated 
on = and = 2. [18-19] 



The persistence of a disturbance of wave vector q at equilibrium is described by 
the normalized time dependent structure factor 
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where the averages on the right hand side of the above mean that is measured in 

equilibrium at time t = then the system is allowed to evolve freely and (j) q is 

measured at time t. Next, we express a At) in terms of the eigenfunctions and 

eigenvalues of the "Hamiltonian" H . Since H commutes with the momentum 
operator, the eigenfunctions of H can be chosen also to be eigenfunctions of the 
momentum operator. Let {y/~ » } be the set of normalized eigenfucntions of H , that 

are also eigenfunctions of the momentum operator with eigenvalue q . The 

corresponding set of eigenvalues of H is (L } . We find that 
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a q (t) = -^- -p— (3) 



It is obvious that if the smallest eigenvalue in the set {X q/3 } , X qo , is really positive, 
oc q (t) cannot decay slower than exponential. Our first goal is to show that to zero 
order in Cl , the volume of the system, X qo is indeed zero at the critical point. 

Consider the set of unnormalized eigenfunctions of the momentum operator with 
momentum q , 

xCq-X X ) = i% i Wg > = q) (4) 



(To simplify we restrict ourselves in the following to sets {^,.}that do not have a 

subset that sums up to zero. This restriction can be avoided in a more detailed 
presentation). The expectation value of H , 



u M(g£ 1 ) _ <X(q-J^-J„)\H\x(q-J l ,-J n )> (5 
<X(qX--J n )\X(qX--J n )> ' 

is an exact upper bound on X qo by the variational theorem (applied to the sector of 
states with momentum q). The expectation value /i can be easily calculated, using 
the identity 



<W G ia - ? # a 9 \Wg ^^Yg I[A_ ? ,[//,A,]]I^g > 



(6) 



The result is 



so that in the limit of the infinite system the sum on the right hand of the above yields 
the exact result. Note, that although the final result may seem to be the result of a 
Gaussian approximation it is in fact the full exact result. The non Gaussian 
corrections vanish when the volume of the system tends to infinity. At the transition 
3^(0) is proportional to £~ 2+T] , so that the term that is zero order in Q. on the right 

hand side of eq.(7) has the form 2_J I,- f I,- = q) .It is straightforward to show that 

if r > 1 such a sum attains its minimal value when 7. =4/ . The value of the sum at 

1 / n 

the minimum is fl (n \q)a n~ {X+e)+n q 2+e ~ n . If 6 is not too negative 1 + 6 -ri is positive 

and fl {n \q) tends to zero as n tends to infinity, for all q . This implies at once that 

X qo is indeed zero. What we have shown is that the (positive) spectrum {X qp } is not 

bound from below by any positive value. We have thus proven that the necessary 
condition for the existence of a stretched exponential decay holds at the point of 
transition. 

We turn now to obtaining a crude estimate of oc q (t) at long times, that is based on 

the simple ideas that have already been presented above. We will assume that the 

system has excitations carrying momentum 7 and energy co t = Bt and that all the 

states of the system can be approximately described as a collection of such 
excitations, so that each state can be tagged by the total number of excitations present 
and the "energy" of a state carrying total momentum q is given approximately by 

X"\q;7 i ,..J„) = £ j a) ei (8) 
i=i 

This assumption is motivated by the form of the trial functions (eq. 5) and the 
corresponding diagonal elements of the Hamiltonian H , (eq. 7). 

The decay function a q (t) is given by 

a q (t) = £}n4 8(£/i -q) \a n (X,-X) I 2 e'^' (9) 

n=l,3 



By the variational theorem it is obvious that z > 1 , so that the minimal value of 
A ( " ) (g;i i ,...,i„) is Bn~ z q z and this value is attained at the point{i ; =qln} . 
Consider the quantity 



K (<?, t) = \ fldlS (£ I -q)\a„ (X ,...,£„) I 2 e^ tf (10) 
1=1 

It is clear that the most dominant contribution to \ (q, t) comes from the vicinity of 
the point {i, = g/n} . Standard asymptotics yields the large t behaviour of A^iqj) 

A„(q,t) = exp[-f(n)-^(n-iyn((O q t)-n- z co il t] . (11) 

In obtaining the above we also required that A n (q,t) = A n {(O q t) . The long time 
behaviour of a q (t) is now determined by obtaining the most dominant A n {(O q t) from 
the equation 

- f\n) - |&i(fiy) + (z " Vn z (O q t = (12) 

Estimation of the function fin) is quite complicated and beyond the scope of this 
paper. It will turn out, however, that much can be said about the decay function 
without knowing the detailed form of fin) . Assume first that fin) is either linear in 

n for large n or that £(jfi = . If this is the case the first term on the left hand 

side of eq. (12) can be ignored and the conclusion is that the dominant contribution 
comes from n = n* , where 

2(z-l) CO J - 

n* = [ y ' — q —Y , (13) 
d £n(O q t 

So that the decay function a q (t) can be written as exp[-F(a> q t)] , where to leading 
order 

F[o)t] a [coty [£ncot] - . (14) 



Apart from the logarithmic correction the stretched exponential decay obtained here is 
exactly the form we found by using entirely different methods, [8,9]. 



Now it is clear that the decay becomes faster when f(n) is increased. Therefore, it is 
not possible that the two point function decays slower than the decay given by eq. 
(14). 



Consider next what happens if f(n) a n a with a > 1 . It is easily verified that in 
such a case the dominant terms on the left hand side of eq.(12) are the first and the last 
and n * giving the main contribution to the sum is 
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n* = [- (O a t] a+ ^ , (15) 

a q 

and 



F((O q t) « [(o q tr +z - x . (16) 



We see that even here the final result is a stretched exponential decay. The stretched 
exponent is larger than before but it is still stretched. Only in the limit we 
recover the unstretched exponential. It is interesting to consider the case where 
indeed the decay is exponential (a = oo)to understand why it cannot be generic. 

Consider the "Hamiltonian" in (eq.2) for the Gaussian model (no 4 term) at its 

transition point. The "excitation energy" 0) q is given by (O q = vqX q and the decay 

function is a q (t) = exp[-co q t] . The reason is that (j) q 1 > is an exact eigenstate of the 

"Hamiltonian". Consequently, the sum in eq. (9) over states with n excitations is 
limited to n = 1 . This sharp cut-off can be interpreted as a = °° . In the non linear 
case (j) q 10 > is far from being an exact eigenstate of the "Hamiltonian" so that it is 

extremely unlikely that a = °° . Furthermore, even if a non- linear model could be 
constructed to yield a = °° , this property would not be stable against small changes in 
the model. Namely a ± °°is the generic case for non- linear systems. In fact, we have 
good reasons to believe that the generic situation is characterized not only by 
a ^ oobut by a < 1 , so that eq. (14) describes the decay correctly. 

We hope to present a detailed estimation of f(n) in the near future. 
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